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Some Consequences of Thermosolutal Convection: 

The Grain Structure of Castings 

G. HANSEN, A. HELLAWELL, S.Z. LU, and R.S. STEUBE 

The essential principles of thermosolutal convection are outlined, and how convection provides a 
transport mechanism between the mushy region of a casting and the open bulk liquid is illustrated. 
The convective flow patterns which develop assist in heat exchange and macroscopic solute segre- 
gation during solidification; they also provide a mechanism for the transport of dendritic fragments 
from the mushy region into the bulk liquid. Surviving fragments become nuclei for equiaxed grains 
and so lead to blocking of the parental columnar, dendritic growth front from which they originated. 
The physical steps in such a sequence are considered and some experimental data are provided to 
support the argument. 


I. INTRODUCTION 

The focus of this article is upon the grain structure of 
castings and whether it is possible to model the columnar- 
equiaxed grain transition in a physically realistic way. The 
argument is that if no heterogeneous substrate materials 
( e.g grain refiners) are deliberately introduced into a melt 
prior to or during pouring, the most probable alternative 
nuclei for equiaxed grains are particles of the base material 
itself. Such crystal fragments may be introduced during the 
pouring operation, especially with low superheats or, at 
longer times, can originate during dendrite ripening within 
the developing columnar mushy region. Without forced 
stirring, the only mechanism for such fragments to “es- 
cape” from the mushy region is by natural convection — 
so-called thermosolutal convection. 

There is nothing new in this argument: it was suggested 
by Rosenhain nearly 70 years ago, [1] and although it has 
been reiterated frequently/ 2 ^ 61 the details of how it might 
operate have been addressed only in a fragmentary way. 
An attempt will be made here to identify the physical steps 
which are implied by this process and to consider how they 
can be linked to provide a predictive sequence for equiaxed 
grain formation and, hence, for the columnar-to-equiaxed 
grain transition. It is in order, first, to outline the charac- 
teristics of thermosolutal convection, which is an essential 
part of the process. 


O. THERMOSOLUTAL CONVECTION 

Three statements, (1) through (3), summarize this phe- 
nomenon. 1 [7 > 8 ’ 91 

1. Natural convection occurs because of density variations 
in a system In a multicomponent system, density dif- 
ferences arise from temperature and compositional var- 
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iations which develop simultaneously. There are there- 
fore two contributions to the buoyancy within a fluid: 
(a) Thermo: 


Pt - Po [• - a (T - r 0 )] [la] 

where p is density, T is temperature, a is the volume co- 
efficient of thermal expansion, which is always positive, 
and p 0 is the density at the melting point, 7 0 . 

(b) Solutal: 


p c = p 0 [l ± /3(C- C 0 )] [lb] 


where the solutal expansion coefficient, /3, can be of either 
sign and p 0 is the density of pure solvent liquid, as de- 
scribed previously. 

Combining these, the density difference between two el- 
ements of liquid, Ap, having temperature and composition 
differences, A T and AC, is then given by 

A p = - p 0 (aAT ± fiAC) [2a] 


or is expressed as a vertical gradient, 


dp 

dz 


dT dC 
p,< “* * f 


[2b] 


Of particular interest will be situations where the signs 
of aAT and (3 AC are opposite, so that it becomes possible 
to have liquid which is both cooler and less dense. In the 
present context, the relevant compositions and temperatures 
are along the liquidus boundary of a phase diagram. 

(2) If an element of liquid translates vertically from one 
position to another in a gravitational field at a velocity 
V (Figure 1), because of an initial density difference 
with respect to surroundings, continued movement will 
depend upon there being retained buoyancy to sustain 
convection and therefore upon how rapidly heat and 
matter exchange laterally between the volume element 
and surroundings. With reference to Figure 1, there are 
two fluxes. 

(a) Thermo: 


J T 




[3a] 


where k and k are the thermal conductivity and diffusivity 
and T and Q are temperature and heat content, respectively 
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gradient is set up around a moving element as it shears 
through the surrounding liquid. If the flow is stream- 
lined, as it generally is in the cases considered here, 
the buoyancy pressure is related to the velocity gradient 
by the dynamic viscosity coefficient, 77: 

dv 

kpgh = 1? — [5] 

dx 

more conveniently normalized as the kinematic viscosity 

coefficient, v = — . 

P 

It is the combination of these three relations that deter- 
mines the scale of convective patterns and rate of flow 
within them. As a general statement, for such convection 
to occur, the diffusivity ratio cannot equal that of the buoy- 
ancy contributions; i.e., 

k , a (dT/dx) 

D l P (dcidx) 1 



Fig. 1 — Schematic for Eqs. [3a] and [b] for a volume element of liquid 
rising at a rate V. 


(k = kiCp , where Cp is the specific heat at constant pres- 
sure); Q and Cp are here defined with respect to volume, 
(b) Solutal: 


J c = D l — [3b] 

dx 

where D L is the liquid diffusion coefficient and C is con- 
centration. 

Of importance here is the ratio of the two diffusivity 
coefficients: 


K 


T 


[4] 


where r is called the Lewis number, having values MO 3 
for metals or MO 2 for aqueous and organic systems. Heat 
redistribution is therefore rapid and the buoyancy is sus- 
tained by a composition difference of much longer duration. 
(3) Fluid flow is opposed by viscous drag and a velocity 


and, in general, r exceeds that ratio by some two orders of 
magnitude. 

In fluid mechanical calculations, it is mathematically 
convenient to describe different situations or problems by 
dimensionless numbers, relating system parameters, veloc- 
ity F(Z,r~‘), and scale (Z,) of events to relevant material 
constants, notably, k, D l , and v(L 2 T ] ). A few of these 
numbers are listed in Table I, with indications of their phys- 
ical meanings and orders of magnitude relevant to the pres- 
ent context. Clearly, the values of those numbers involving 
system parameters are particularly sensitive to the choice 
of that dimension, L, describing the scale of convective 
flow, e.g., height or width of the system or of the local flow 
pattern. Since, in the present context we are concerned with 
local flow plumes of the order of 10 _3 m wide (but much 
longer), the numbers quoted are for that scale. We note in 
particular that the Reynolds number, even for metallic sys- 
tems, is below the limit (~*10 3 ) at which streamlined flow 
breaks down into turbulence. It will also be apparent that 
Peclet numbers, Pe r or Pe^ are products of the Reynolds 
number and the Prandtl or Schmidt number. Keeping the 
physical significance of these numbers in perspective is cru- 
cial to understand their implications. Thus, referring to the 
Peclet numbers, a value of one would indicate that heat or 
mass transfer took place as rapidly as convective flow, so 
that the fluid would remain quiescent. In the present con- 
text, the two Peclet numbers indicate that fluid flow typi- 
cally takes place a little faster than allows thermal 
equilibrium to be attained, but very much faster than com- 
positional (solutal) adjustment, cf Eqs. [3] and [4]. The rel- 
atively big difference between Prandtl numbers for metals 
vs aqueous and organic systems reflects their higher thermal 
conductivity, but as it transpires, this does not significantly 
affect the scale of convective patterns, although the flow 
rates are much rapid. 


HI. PHENOMENOLOGICAL CONSEQUENCES 

Figure 2 shows three schematic casting configurations 
and the form of thermosolutal flow patterns which are to 
be expected or are observed. The left-hand sequence, I 
through III (a), corresponds to dpldz negative, Le., solute 
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Table I. To Identify Relevant Dimensionless Numbers 


Name 

Ratio 

Physical Meaning 

Order of Magnitude 

Prandtl 

Schmidt 

a = v!k 
SC = v/D 

velocity vs temperature field 
velocity vs composition field 

— 10 2 (metals) to —10 (aqueous, organic) 
— 10 2 (metals) to — 10 3 (aqueous, organic) 

Reynolds 

£ 

ii 

* 15 

convective vs diffusive momentum transfer 

—5 X 10 2 (metals to < 10 2 (aqueous organic) 

Peclet 

no -o 

r& o 

o, -H 

II II 

15 

convective vs thermal transfer 
convective vs solutal transfer 

-10 

— 10 4 (metals) to — 10 3 (aqueous, organic) 


sufficiently less dense to outweigh the thermal buoyancy 
(cooler liquid is less dense), and the right-hand sequence, 
I through III (b), is the reverse, with more dense solute. 
These would correspond to lead base alloys with tin and 
vice versa . The shaded regions correspond to columnar 
mushy regions, growing into positive temperature gradients 
with the gravitational vector vertically downward. The ar- 
rows, K, indicate the directions of growth. Cases 1(a) and 
(b), with growth vertically upward, have the supernatant 
liquid stabilized against convection by the positive temper- 
ature gradient, but there is density inversion in case 1(a) 
because the cooler interdendritic liquid is less dense; this 
will be of central interest hereafter. In case 1(b), the entire 
liquid is quiescent; dpidz is negative. In cases 11(a) and (b), 
the open liquid will be mixed by thermally driven convec- 
tion; in case 11(a), the less dense interdendritic liquid will 
not contribute to convection, but in case 11(b), the solutal 
contribution will tend to cause downward flow from the 
growth front, interacting with thermal convection. In cases 
111(a) and (b), with solidification from the sides, solutal 
buoyancy will cause (a) upward or (b) downward convec- 
tion at and behind the growth fronts, and these will also 
interact with thermal convection which is not here repre- 
sented. Macroscopic solute segregation in these two cases 
will typically reduce the freezing point at the (a) top or (b) 
bottom of the ingots, so that the upper and lower parts of 
the growth fronts present lopsided versions of 1(a) or 11(b). 

We will concentrate on the configuration of 1(a), because 
it allows the thermal and solutal contributions to be most 
easily separated and is thus the most analytically tractable 
of these alternatives. References to the other configurations 
will be made occasionally, but in any case, extensions to 
these can easily be anticipated in a qualitative way. 

Referring then to case 1(a) (Figure 3), from left to right: 
(a) alloy of composition C 0 in a phase diagram with k 0 < 

1 and negative liquidus slope, m L , grows upward with liquid 
concentration rising to C f at dendrite tips; (b) there is a 
concentration ripple across the front of the wave length de- 
termined by the primary dendritic spacing; (c) concentra- 
tion profiles along the vertical axes A and B then appear, 
with interdendritic liquid, trace B, rising according to the 
reciprocal liquidus slope and vertical temperature gradient, 
i.e., dCIdz = dddT • dTidz\ and (d) the corresponding ver- 
tical density profiles showing, with /3 dCIdz > — a dTidz 
(i.e., Eq. [2b]), a small dip at dendrite tips, corresponding 
to C„ while the interdendritic liquid shows density inver- 
sion with a positive value of dpidz. (If the solutal coefficient 
/3 = 0, the density profile would continue downward ac- 
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(a) (b) 

Fig. 2 — Alternative casting geometries, I through III, with schematic 
thermosolutal convection patterns, solid lines within open liquid, broken 
lines; interdendritic: shaded regions correspond to permeable mushy 
regions: (a) for cases where the solute is less dense than the solvent, 
fidCldz > —adTldz (cooler liquid is the less dense) and (6) for the reverse 
situation where the thermosolutal buoyancy terms have the same sign 
(solute is the more dense). The arrows, K, indicate the growth directions: 
I — upward, II — downward, and III— from the sides, horizontally. The 
gravitational vector, g, is vertically downward in all cases. 


cording to the temperature gradient.) With such density in- 
version, the system is unstable and thermosolutal 
convection is expected. The critical region for the pertur- 
bations) leading to such convection is close to the growth 
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dp/dz=-pObdT/dz 


Co 



(a) (b) (c) (d) 

Fig. 3 — (c) Composition and (</) density profiles through ( b ) an upward growing dendritic front (a) Trace A along a dendrite axis and trace B, interdendritic, 
with reference to a phase diagram and alloy composition C. 


front where the density gradient changes sign, and experi- 
mentation confirms that this is so. [1 2,131 Perturbations at 
lower positions within the mushy region are improbable. 
The consequences of this configuration are well docu- 
mented for transparent system s {1CM51 and are briefly as fol- 
lows: very early (CIO 2 s) in the growth of the dendritic 
front, short-range convective “fingers” develop, rising up 
to **10 mm into the supernatant liquid with concomitant 
downward flow to create a recirculating, mixing zone, 
which is slightly solute enriched, approximately 1 to 2 pet. 
At somewhat longer times (> 10 3 s), finger convection cou- 
ples with the interdendritic liquid and locally reinforces to 
produce longer range (>100 mm) solute-rich plumes. As 
these develop, the shorter range patterns are largely re- 
placed and, simultaneously with plume formation, channels 
develop within the dendritic array below. These features 
can be observed directly in a transparent aqueous or organic 
system and Figures 4 (a) through (c) are examples from the 
aqueous ammonium chloride water system, in which water 
is the less dense solute. Figure 4(d) is an example of such 
c hann el formation in a metallic ingot, where accumulation 
of dendritic debris within channels gives rise to character- 
istic polycrystalline freckles which delineate the original 
events. [151 In Figure 4(a) is shown the short-range, oscilla- 
tory finger convection patterns, while Figure 4(b) shows the 
longer range, nearly steady state plumes, which rise from 
channel mouths, such as that shown in Figure 4(c). The 
dimensions of channels and dendrites are very similar in 
aqueous and metallic systems, although the flow rates in 
the latter are undoubtedly much fasterJ 151 

Measurement and analysis of channel-plume motion in 
transparent systems allows extrapolation to opaque metallic 
systems, based on postmortem examination of the fully so- 
lidified samples^ 14 - 151 Plume widths are all about 10~ 3 m 
(some three to five primary dendrite spacings), and plumes 
are enriched in solute by 1 to 2 pet compared with the bulk 
liquid; as liquid leaves the dendritic front, it is also a little 


cooler (<5 K) than the surroundings but rapidly warms up 
to the ambient temperature (c/ thermal Peclet number). 
Flow rates in the plumes have been measured in the trans- 
parent materials and analysis of these, extrapolated to met- 
als, indicates mean velocities, organic:aqueous:metallic, of 
about 10" 3 ms -1 : 10 -2 ms -1 : 10" 1 ms~\ the last being very 
rapid but still streamlined. At a given time, some 2 to 3 pet 
of the liquid flows upward rapidly in such plumes, while 
the balance is drawn down into the upper part (<10 mm) 
of the mushy region to be continuously entrained in the 
plumes. This means, as shown in Figure 4(c), that solvent- 
rich, warmer liquid is drawn through the upper part of the 
dendritic array. The finger convection also involves some 
short-range recirculation between dendrites, but it does not 
penetrate the mushy region to greater depths than the pri- 
mary spacing (<10~ 3 m) and there are no channels asso- 
ciated with it It may be noted that in a configuration with 
horizontal growth, as in a billet casting. Figure 2, III, the 
finger zone can never be fully established but is mostly 
swept away by upward or downward flow across the growth 
front. Plume flow is similarly obscured at such a front, but 
this does not mean that exchange between the mushy region 
and open liquid is absent. 

In all cases, the convection currents carry entrained den- 
dritic fragments within them and thus provide particle trans- 
port over shorter or longer ranges into the bulk liquid. This 
will be our principal focus. 

It should be added here that there are two other conse- 
quences of this convection, namely, macroscopic thermal 
and solutal mixing within the open bulk liquid. These are 
particularly obvious in the vertical configuration of Figure 
2, 1(a). The re-entrainment means that there is accelerated 
heat exchange between the open liquid and the mushy re- 
gion, causing the former to cool more rapidly than it would 
by diffusion alone. Simultaneously, liquid mixing causes 
macrosegregation and progressive depression of the freez- 
ing point, so that the schematic concentration profile of Fig- 
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3 cm 


(d) 

Fig. 4— (a) Side view of short range finger convection above a dendritic growth front in NH 4 CI-H,0; (6) long-range plume convection from channel 
mouths at a similar growth front. The optical contrast in both cases arises from composition-dependent variations in refractive index, (c) Detail of a channel 
mouth from which a convection plume, (b), is flowing, showing scale of channel and dendrites: and (</), surface freckles on the outside of a superalloy 
ingot (courtesy of Giamei, United Technologies Research Center, Dayton, OH). Scales as shown. All relate to case 1(a) of Fig. 2. 

ure 3(c) is not actually steady state, i.e., C L # C 0 . During figuration. With longer range mixing of the liquid, the sol- 

the regime of finger convection, growth takes place into a ute redistribution can be described by an effective 

solute-enriched zone, comparable with a zone-melting con- distribution coefficient, k cff > where < k etT < 1. Without 
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Table II. Model Sequence for Grain Structure (Events/ 
Volume/Time) 


(1) 

Side arm generation and ripening 
(>10* cm -3 at >10 2 cnrV 1 ) 



(2) 

Side arm separation with and without fluid flow 
(fraction >10 ') 


i 

(3) 

Fragment escape by convection/stirring 
(>10 cnrV) 


i 

(4) 

Fragment survival 

(melting/life span/population density) 


i 

(5) 

Growth and sedimentation 
(modified stokes) 


i 

(6) 

Columnar blocking and grain size 
(columnar/equiaxed) 


such liquid mixing, steady state prevails and k eff = 1 ; with 
mixing, k eff may be as low as 0.9 and the mean mushy zone 
composition is adjusted by a larger solid fraction than the 
lever rule would predict. This has been analyzed else- 
where. 1131 

With these essential features in mind, we proceed to the 
question of modeling how the grain structure of a casting 
can develop when no relevant heterogeneous substrates are 
present in the liquid. This is the situation if deliberate ad- 
ditions have not been made to the melt and assumes that 
dendrite fragments of the base material must be better “nu- 
clei” for equiaxed grains than any other accidental “dirt” 
in the system, such as oxide or slag particles. 

IV. A PHYSICAL MODEL 

In Table II is a model to describe the development of an 
equiaxed grain structure from intrinsic nucleation sites in- 
volving sequential steps (1) through (6), as follows: 

(1) Side arms of dendrites are continuously generated and 
coarsen by ripening behind the columnar growth front; 

(2) A fraction of these side arms become detached from 
primary dendrite stalks; 

(3) A further fraction escape from the mushy zone into the 
bulk liquid; 

(4) Some of these survive, subject to time, local tempera- 
ture, and composition; 

(5) Survivors grow in supercooled liquid and sink as they 
do so, 

(6) Survivors eventually block the advance of the parental 
growth front from which they originated. 

The details could be further subdivided according to 
choice, but this is a closed system (any casting, even a 



Fig. 5— NH 4 C1 dendrites growing horizontally on a gradient stage: (a) 
under steady-state conditions, with little, if any, secondary arm 
detachment; and (b) after a small temperature rise, when nearly all side 
arms become detached. 

continuous casting, constitutes a closed system) and the se- 
quence describes a physically realistic and observed situa- 
tion. There are no nucleation kinetics in the classical sense, 
apart from those of a chilled layer on the mold walls, from 
which the columnar, dendritic front evolves. We consider 
these steps separately, indicating order of magnitude esti- 
mates where or when available, referring mostly to obser- 
vations of the transparent system NH 4 C1-H 2 0. 

A. Dendrite Structure 

Figure 5 shows NH 4 C1 dendrites growing horizontally in 
a thin cell (—125 jjl m) on a temperature gradient stage. 
There could be little or no convection in such a cell. The 
primary dendrite spacing is comparable with, or a little 
larger than, that observed in a casting of the same nominal 
composition. Many similar pictures have been shown in the 
literature/ 416171 and we will make only a few observations. 
The density of side arms is very high and in a “typical” 
mushy region equates with some 10 6 cm -3 within this re- 
gion. The number of potential crystal fragments is enor- 
mous, even disregarding tertiary arms. With a growth rate 
of — 10~ 3 cm s~ l , side arms are continuously generated at 
a rate around 10 3 cm -2 s _l across the dendritic growth front. 
Side arms have initial, intrinsic, or nascent spacings. A?, 
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which are in a simple ratio to the radii of curvature, r, of 
the primary dendrite tips and hence can be related to the 
imposed growth rate, composition, temperature gradient, 
and primary spacings, as analyses have shown/ 18 * 231 In this 
material, over an order of magnitude in the growth rate, 
this ratio is approximately A ° 2 :r = 4:1. We take this as a 
reflection of the oscillatory growth of dendntes, growing at 
the condition of “marginal stability,” and believe that side 
arms grow from the initial perturbations so created at the 
primary tips. Discussion of the geometry, wave length, and 
frequency of dendrite tip perturbations is beyond the im- 
mediate scope of this article, but we note that side arm 
spacings are always in a simple ratio to the tip curva- 
ture, 119 * 201 whether in single-component systems (heat flow 
only) or with solute present. Ripening of the intrinsic spac- 
ing begins almost at once and proceeds according to a re- 
lationship of the form A 2 - A ° 2 = Kr , where t = time and 
n lies between 1/3 and 1/2 depending on the solid volume 
fraction, as described elsewhere in this series/ 211 This step 
is sufficiently well documented to be predictive. 

B. Side Arm Detachment 

It may be seen from Figure 5(a) that very few side arms 
have separated from the primary spine, although the junc- 
tions between them have necked down to small radii of 
curvature. This seems to be the situation if the temperature 
gradient is constant and the sample is cooling continuously. 
However, a small temperature rise (as yet undetermined) 
causes nearly all the side arms to detach abruptly (Figure 
5(b)), as Sato et al. [22] also demonstrated with succinonitrile. 
In practice, such a temperature rise can occur if the mold 
temperature rises, as in a sand casting, or, of particular rel- 
evance, when the growth rate decreases, as it does contin- 
uously in the columnar region of a casting. With 
deceleration, there is an immediate increase in the dendrite 
tip radii, but a much slower increase in the primary spac- 
ings (e.g., Figure 7 of Reference 4). Consequently, the in- 
terdendritic solute concentration is raised above that for 
growth at the former growth rate and this causes side arms 
to melt off at their narrow roots. With such transient growth 
conditions, dendrites are always a little closer together then 
they would be under steady-state conditions/ 401 

Better temperature measurement and analysis of this sit- 
uation are needed if there are to be predictive statements. 
We note that ripening, detachment, and re-entrainment take 
place continuously and concurrently and that the permea- 
bility of the dendritic array is also changing with position 
behind the growth front. 

C. Fragment Escape 

Reference to the transport of dendrite fragments by ther- 
mosolutal fingers and plumes has already been made. In the 
vertical configuration of Figure 2, 1(a), finger convection is 
of short range and involves only local recirculation. An 
estimate from video records would be around 10 particles 
cm -2 s _1 so ejected, having diameters between 10 to 40 
pum. In the vertical configuration, these may not be of great 
importance for future equiaxed grain formation, but in a 
casting with side chill, such currents are swept away into 
the open liquid, as already noted. Within plumes, fragments 
are more numerous and are traveling faster, but with the 


spatial separation of channels of ^5 mm, this again equates 
with a mean rate of escape around 10 to 30 fragments cm' 2 
s~ l , against some 500 primary dendrites cm -2 . With the 
regeneration of side arms, around 200 cm' 2 s ', previously 
noted, this must mean that about 10 pet of the potentially 
available fragments leave the mushy region and therefore 
that the fraction detached must be greater than this, but we 
do not yet know by how much. For the present, in the 
aqueous system, it will be safe to assume a particle escape 
density of >10 cm" 2 s _I , with a lower size estimate limited 
by resolution to «*10 pm or, on average, —20-yu.m diame- 
ter. There can be no estimates for a metallic system, as yet. 
Flow rates in metals should be an order of magnitude 
greater within plumes of similar diameters and spacings/ 151 
We are limited here by uncertainty about the detachment 
kinetics. 

D. Fragment Survival 

Three areas, (a) through (c), may be identified: 

(a) To treat the problem of survival requires a detailed ther- 
mal record with respect to time and position, in this 
case, as vertical temperature profiles, against which are 
required the corresponding particle population densi- 
ties. Schematically (Figure 6(a)) the bulk liquid is cool- 
ing by heat extraction from the chilled base, assisted 
by convective mixing and heat exchange with the 
mushy region during re -entrainment. The temperature 
profiles are hinged about the growth front at T h so that 
predictive analysis of the interfacial undercooling, A T h 
is an essential part of this problem (^References 18 
and 23). Temperature gradients in the liquid rapidly de- 
crease with time, especially in metals and with convec- 
tive mixing, so that in a very short time interval (10 to 
10 2 s), there is a transition from most of the liquid being 
slightly superheated to its being almost uniformly su- 
percooled. Above the freezing point, fragments melt, 
but as the profile falls below T Ly all survivors start to 
grow, recalescence occurs, and the liquid undercooling 
becomes small and uniform. With knowledge of the 
material properties (conductivities, specific, and latent 
heats), the liquid composition and, with a model for 
undercooling at the growth front, the temperature pro- 
files can be computed with readily available thermal 
codes. This is a simple one-dimensional analysis which, 
however, must also incorporate thermal mixing in the 
open liquid by convection and the attendant heat 
exchange during re -entrainment. Experimentally, the 
temperatures can also be followed by a series of ther- 
mocouples or, better, can be scanned by one thermo- 
couple, moving vertically up and down through the 
melt. Figure 6(b) is an example of such a scanning se- 
quence for an aqueous NH 4 C1 casting, poured into a 
mold with an initial superheat of *»30 K above the li- 
quidus temperature. The positions and temperatures of 
the growth front are indicated on these temperature pro- 
files and reflect the depression of the growth front from 
changes in bulk liquid composition (the finger zone 
dominates up to a time of about 30 minutes, after which 
macrosegregation occurs) and from the front under- 
cooling, AT, 3 to 4 K. 

(b) Until the bulk liquid falls below the freezing point, 
ejected particles melt and have a limited life span. The 
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Temperature Profiles 


— ■ — 2 min. 


8 min. 


18 min. 


27 min. 


-* — 30 min. — A — 44 min. — * — 49 min. — 0 — 60 min. 



Ti To 

Pig (y — (q) Schematic vertical temperature profiles, G ; through G 3 , at successively later times, above a dendritic front, growing at an undercooling, AT), 
for velocity, K, and ( b ) actual temperature profiles recorded by a single thermocouple, scanning vertically through an NH<C1-H 2 0 casting (positions of the 
growth front are shown by arrows). 


kinetics of melting depend upon the ambient tempera- 
ture, composition, and rate of movement of the particle 
through the liquid. In the present context, the total mass 
of solid is very small by comparison with the volume 
of liquid, so that the latent heat required for melting 
and die composition changes resulting therefrom are 
negligible. The problem may therefore be described as 
one of “unconstrained*’ melting. The rate of melting 
will clearly be a sensitive function of movement with 
respect to the surroundings, whether by long-range dif- 
fusion (of heat or matter) or through a boundary layer 
of limited thickness, but initially, it is necessary to sep- 
arate the contributions from heat and matter transport. 
It will be assumed here that the times involved are 
small and that solid-state diffusion is not important. Ex- 
perimental observations of ammonium chloride parti- 
cles in aqueous castings indicate that in a thermally 
hostile environment, fragments of 10- to 20-/xm di- 
ameter rarely survive longer than ^10 seconds. The 
term “hostile,” of course, needs qualification, because 
in the finger regime, there is continuous mixing of 


cooler, water-rich liquid with warmer salt-rich liquid, 
and plume temperatures are initially below those of the 
bulk liquid, so the precise environment is difficult to 
specify. However, reference to the schematic phase di- 
agram of Figure 7 will show that there must be a tran- 
sition in melting kinetics from a regime controlled by 
matter transport to one dominated by heat flow. 
Consider (Figure 7), a bath composition of C 0 on a phase 
diagram of given liquidus slope (solid line of slope m L = 
— A77AC). In the first place, consider a system such as 
NH 4 C1-H 2 0 having zero solid solubility, = 0. A particle, 
initially at equilibrium at T 0 , is placed in a large volume of 
liquid, C 0 , at a superheat of A T hy where the equilibrium 
solute concentration is lower, by AC. Melting is then at a 
rate V m , given by 


V m (C, - C s ) 


DtAC 

Ar 


[7a] 


where C, is the liquid interfacial composition at T 0 + A T h , 
C s is the solid composition = 0, D L is the liquid diffusion 
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Liquidus 


Solidus 

Thermo/So lutal Transition 
Temperature of Melting 

Fl 8; 7 ~ T ^ sltlon from solutally controlled to thermally controlled melting at a temperature, T ML> for a case of zero solid solubility, ^ = 0 and finite 
solid solubility, ^ - 0.25. The solid line is the liquidus, and the bath composition is C 0 , with freezing point 7;0. At temperatures above tJ up to 7V.. 
melting is controlled by solute diffusion and is slow; above T ML , melting depends on heat flow and is rapid. The melting rate, V , at any temperature up 
to is proportional to AC « C 0 - C ; . " K p 


coefficient, and A r is some radial diffusion distance in the 
liquid, e.g. y a boundary layer, ~5. The rate is limited by 
solute/solvent diffusion and rises with the temperature as 
C { decreases and AC increases, until T ML . At and above the 
melting point of the pure solvent, melting becomes heat- 
flow dependent and reduces to a single component problem; 
composition adjustment then follows after or while the 
particle melts, which it does increasingly rapidly (cf r ~ 
k/DJ as 


V m AH = 


*IL I TJ 

A r' 


[7b] 


where T is a bath temperature above the melting point, A H 
the latent heat, k the thermal conductivity, and Ari a larger 
characteristic distance for heat flow. In a more general case, 
with some solid solubility, e.g = 0.25 (Figure 7), the 
same argument applies, but now the transition temperature 
from matter- to heat-controlled melting applies to a bath 
temperature for liquid C 0 at the corresponding solid melting 
composition at k 0 C 0 . The locus of such temperatures then 
follows a line 


= T m ~ m L KC [ 8 ] 

It seems appropriate to call this a thermosolutal transition 
temperature for melting. Crude experiments with succinon- 
itrile-water confirm changes in melting rates above T mt 
(-Tm) by a factor of more than 30 times within less than 
5 K above the melting point. r241 For the NH 4 C1-H 2 0 system. 


this transition is irrelevant, because the aqueous solutions 
are far below the solid melting point, but in foundry prac- 
tice, in many metal systems, the pouring temperature may 
often exceed this melting transition temperature (for ex- 
ample, consider alloys such as Fe-0.5 wt pet or Al-4 wt pet 
Cu). This will mean that any metal fragments which are 
momentarily produced, as by local chill at mold surfaces, 
will have a very short lifetime indeed, unless the superheat 
is almost impracticably low. Reverting to the experimental 
aqueous system under consideration, approximate melting 
rates and hence life spans of the same order as are observed 
can be calculated for solute/solvent diffusion by taking a 
boundary layer thickness equal to the particle radius. These 
allow some estimates of expected particle densities to be 
made for given local superheats. 

(c) Particle densities prevailing at a given ambient temper- 
ature must be a result of supply and survival, being a 
standard demographic problem of “birth rate” vs 
“death rate,” where here, the former is the number of 
particles supplied by convection, ft and the latter is the 
reciprocal of the melting life span, <f> s _l . The number 
of particles surviving in a given volume after a time, /, 
is then given by 

dft 

— - N - N t <f>, and with integration, [9] 

= ( -j H - exp (-<£f)]) 

<P 
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Pig 8 — Calculated particle population densities, N n after times, f, for a 
given rate of particle arrival or birth rate, jV, and of particle radius, r q, at 
melting or death rates, <f>> corresponding to various melt superheats, A 7^, 
as in Fig. 7. 

Experimentally, we do not yet have acceptable meas- 
urements of particle densities or of the quantities, ft and 
<j>, although an initial value for ft can be taken from 
numbers observed in fingers or plumes and estimates 
of <f> are at least of the correct order. Thus, Figure 8 
shows how particle densities should converge with time 
for assumed initial particle size (r 0 = 10 /xm) and an 
observed supply rate of the order of 10 7 m~ 3 s _1 at 
superheats from 1 to 10 K, Convergence is more rapid 
as the melting rate or death rate increases, but the val- 
ues are considered in general to be too high because 
many particles melt as they are transported through the 
hotter liquid before reaching some given volume ele- 
ment, as previously mentioned; i.e.> ft is probably too 
high. Clearly, the situation of particular interest is with 
the small critical temperature range as the superheat 
falls below — 1 K, when nearly all ejected dendritic 
fragments begin to survive. Meaningful measurements 
of particle densities require very close correlation with 
accurate temperature measurements and local cooling 
rates, and the important transition range is narrow and 
almost abmpt. 

E. Growth and Sedimentation 

Below the freezing point, the situation changes abruptly, 
and from the foregoing (Figure 6(a)), this applies to much 
of the bulk liquid above the growth front. Initially spherical 
particles start to grow dendritically as soon as the under- 
cooling exceeds that for curvature, A T r (with r = 10 /im, 
a *** 0.1 J*m -2 and the entropy of fusion, A S ** 4.10 8 J*m” 3 
K" 1 , and &T r = 2crlrAS < 10“ 4 K). Initially, this will be 
unconstrained growth and will lead to recalescence within 
a short time, so that the bulk liquid rapidly becomes spa- 
tially isothermal. Simple calculations of unconstrained 
growth with solute rejection and with a cooling rate of 
«*10 _1 K # s -1 , from a given number of intrinsic nuclei, in- 
dicate that recalescence will occur within <10 seconds. Re- 
calescence is a consequence of growth and evolution of 
latent heat, so that within quite wide ranges of nucleation 
densities, the rate of recalescence does not actually differ 
very much, perhaps contrary to perceived wisdom. 1291 At 
this point, it is recognized that the modeling procedure be- 
gins to follow that which is generally used to fit with ther- 


mal analysis, cooling curve data [25 281 and becomes 
“standard.” As noted previously, relatively close correla- 
tion with cooling curves is possible with a relatively wide 
range of nucleation sites. 

In this range (of time and temperature), before impinge- 
ment, the radial dendritic growth rate, follows: 

V r « AT” [10] 

where n « 2, and experimental data and analysis are avail- 
able for this situation 116,17,201 in single- and two-component 
systems. 

Since the solid is invariably more dense than the liquid 
(by some 30 pet in NH 4 C1-H 2 0), the growing equiaxed 
crystals sink according to Stokes’ law, modified for an en- 
velope of radius r, within which the solid fraction, f 5 , is 
relatively low (<5 pet). Terminal sedimentation velocities, 
Vj, are reached rapidly in the density ranges involved here, 
although the particles accelerate as they grow. Also, within 
the envelope which effectively encloses the dendritic crys- 
tal, the interdendritic liquid is slightly enriched in solute 
and is therefore marginally less dense than that of the sur- 
rounding liquid. However, although this fraction is large, 
the density difference, Ap r/L , is very small and the buoy- 
ancy contribution is negligible. The terminal velocity, 
may therefore be written 

V T = otg [/jA p s/L + (1 - fs) A p LV J r 2 /rf [11] 

where a < 2/9 , g is the gravitational acceleration,^ is the 
solid fraction within the grain, Ap^ L is the solidrbulk liquid 
density difference, and A p LVL is the difference between the 
liquids inside and outside the envelope. The value of a is 
2/9 for a sphere and the dendritic crystal outlines are ap- 
proximately octahedral, but small shape changes make sur- 
prisingly little difference to Stokes’ law, 1301 and in any case, 
crystals often rotate as they fall. Of more importance is the 
interaction between falling crystals, and as the separation 
between them decreases below about 10 diameters, the ve- 
locities decrease because of an effective upward displace- 
ment of the medium. 1311 It follows that as the particle 
density and size increase, the regime must change sharply 
as impingement begins, from one of falling particles 
through relatively stationary liquid to one of fluid flow 
through a stationary dendritic mesh. 137,381 Analysis of sedi- 
mentation rates in the NH 4 C1-H 2 0 system gives good agree- 
ment with observations with a *** 10 -1 and f ** 0.02 and 
has been scaled with some success to estimate rates and 
grain sizes in a continuous steel casting. 1311 Figure 9 sum- 
marizes this situation schematically, bringing the cycle of 
operations to a conclusion with the blocking of the colum- 
nar front from which fragment nuclei originated. Figure 10 
is an example of this situation in the NH 4 C1-H 2 0 system 
prior to the transition. 

F. Columnar-Equiaxed Transition 

The final step concerns a criterion for how many falling 
equiaxed crystals, and of what size(s), are necessary to 
block an upward growing, columnar dendritic front. Refer- 
ring to Figure 10, a crude estimate of this necessary block- 
ing density might be made as follows. 

The columnar front is growing with a velocity, F c , and 
with a primary dendritic spacing, A t , which is typically the 


578 — VOLUME 27 A, MARCH 1996 


METALLURGICAL AND MATERIALS TRANSACTIONS A 



[ 12 ] 


Sedimentation-Columnar/Equiaxed Transition 



Radial Growth: 

v r <~Ar 

Diameter: 

docjvdt 

o 

Terminal Falling Rate: 

Time to Fall Through H ; : 

/ oo HJV t 

Arrival Rate Needed to Block: 

N a oo VJdt ( Area ‘Time 1 ) 


Fig. 9 — Schematic illustrating the growth of dendritic equiaxed crystals 
of diameters, d y falling at terminal velocities, V ^ toward the dendritic front 
which is growing upward at a velocity 



100 [im 


Fig. 10 — Equiaxed crystals of NH*C1 falling onto the dendritic growth 
front in an aqueous casting. 


same or less than the diameter, d , of falling crystals. The 
basic requirement to avoid blocking is that the columnar 
front be able to grow through a broken layer of obstacles 
in less time than it takes another equivalent layer to arrive, 
so that by dendrite branching and the development of ter- 
tiary arms, as indicated in Figure 9, the primary structure 
is largely retained. A simple criterion might be that the 
critical coverage of the surface should be ^O.S in this time 
interval. Whatever the exact real fraction, the rate of arrival 
of obstacles per unit area, N A , of diameter, d , will then be 
given by 



(where K = 1/4 for 0.5 fraction). There are some obvious 
simplications here, notably, that all the falling crystals will 
not be the same size (coming from different heights) and 
that they also grow after they have arrived on the front 
(mostly in the vertical direction). A more comprehensive 
and sophisticated analysis* 331 is to treat the arriving crystals 
in terms of nucleation and growth according to an Avrami 
formulation for extended volumes, and this actually does 
indicate a necessary arrival rate to effect *=»(). 5 coverage. 
For all these quantities, dimensions and rates can be meas- 
ured and analyzed with some confidence; two points may 
be made. First, the blocking is obviously more effective by 
larger crystals, so that those originating only a little above 
the front are not as effective as those which have fallen 
further and had more time to grow (and also arrive faster). 
Therefore, there is a “most probable” level of origin for 
those equiaxed nuclei which are the most likely to block 
the front. Application of such analysis has been made for 
the model system and can be extrapolated with some con- 
fidence.* 32,331 Second, unless there is a sudden “blizzard” 
of equiaxed crystals (which is often the case, however), the 
columnar-to-equiaxed transition is not abrupt but rather 
gradual, as stray equiaxed grains become incorporated into 
the general framework of the columnar array. Again, this 
is rather well recognized and observed in practice. 


V. SUMMARY AND CONCLUSIONS 

In the foregoing, we have emphasized the roll of ther- 
mosolutal convection as a mechanism to transport ciystal- 
line fragments within a casting, in particular, from within 
the mushy region of a columnar dendritic array into open, 
bulk liquid. The argument pursued throughout has been that 
such fragments are the most probable nuclei for equiaxed 
grains if no other substrates are deliberately introduced into 
the system. 

Present schemes for the modeling of the grain structure 
of castings follow, with some sophistication, that of Old- 
field,* 251 in which the question of origin of nuclei is avoided 
by use of expressions based on classical nucleation kinetics 
upon unspecified sites of variable potency. The modeling 
procedures which are then used, however elegant and elab- 
orate they may be,* 28,34 ’ 351 must be classified as curve-fitting 
exercises to cooling curve data. This is not to detract from 
the demonstrable achievements of these models, but the 
present assertion is that sites can be identified in a physi- 
cally realistic way and the separate steps to describe this 
alternative approach have been considered. This is not a 
simple undertaking and the necessary confirmation from ex- 
perimental measurements is lacking for some parts of the 
operation, although the physical processes can be identified, 
albeit in a somewhat qualitative way. 

Referring to the possible modeling sequence of Table II, 
the origin of crystal fragments is clearly identifiable (1), as 
is their transport (3), but the actual detachment kinetics of 
side arms from primary dendrite stalks (2) have yet to be 
clarified and experimentally confirmed. Likewise, working 
back from the end, reasonable physical understanding of 
the columnar blocking (6) and of the supply of equiaxed 
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“obstacles” (5) is available, supported by acceptable ob- 
servations and measurements. But, at present, the central 
connection, in terms of particle survival and particulate 
density (4), seems to be less tractable, although the kinetics 
of melting appear to be relatively straightforward, if poorly 
documented. 

Inasmuch as a casting is a closed system, it would, in 
principle, be possible to write a computer program for 
stages (3) through (4), the number of particles escaping to 
that number needed to effect blocking, and to simply for- 
malize the intermediate steps to yield a numerically com- 
patible connection. Wang and Beckermann [361 have outlined 
such a scheme to do this, but there is the danger that the 
physics of each individual, concurrent step will become ob- 
scure and that the model will then become no more phys- 
ically predictive than those currently employed. This could 
easily become a more cumbersome and invidious exercise, 
even omitting the dendrite ripening and detachment proc- 
esses and, as such, would offer no significant improvement 
over current models. However, we conclude that despite 
some current difficulties, solutions to the overall modeling 
sequences must be attainable, if only because the process 
is physically correct. It is relevant to note that in practice, 
the ranges of temperature which allow a sound casting to 
be made, combined with the physical properties of mate- 
rials, are quite narrow, and in practice, without the addition 
of grain-refining substrates, the observed range of equiaxed 
grain sizes is limited to not much more than one order of 
magnitude, from around a hundred microns to a few mil- 
limeters. There is not much latitude here for meaningful 
results, although in terms of nucleation density, this size 
range does correspond to more than three orders of mag- 
nitude and is not as narrow as the grain size itself might 
imply. 

Some comments are in order about the value of using 
transparent aqueous or organic analogues for metallic be- 
havior and about extension from the vertical unidirectional 
casting geometry to other more general three-dimensional 
configurations. 

There are two major inter-related differences between 
aqueous and metallic systems which require careful scaling 
because of the relative thermal conductivities and temper- 
ature ranges typically involved. The first of these is re- 
flected in the Prandtl numbers (Table I) which, for metals, 
are some two or more orders of magnitude smaller than for 
the transparent materials. As previously remarked, this does 
not appear to alter the scale of convective patterns, but it 
profoundly affects the flow rates within them. Simultane- 
ously, however, the higher thermal conductivity also allows 
heat flow to keep up with the rapid transport over similar 
distances (but shorter times) to those in the other systems, 
so that estimates of the distances within which convective 
plumes approach thermal equilibrium with surrounding liq- 
uid are quite similar, around 10 mm. [15) The scales remain 
similar but the motion is dramatically accelerated, and there 
is some doubt as to whether a well-defined finger regime 
has time to develop in metals. In a base-chilled, lead alloy 
ingot, vertical segregation channels can be observed within 
a few millimeters of the chilled face* 391 and there is a very 
brief incubation period compared with aqueous systems. 
We conclude that in this respect, the transparent materials 
are valid analogues but in “slow motion.” 


Brief references were made to extrapolation from the ver- 
tical to other configurations, where thermosolutal and/or 
purely thermal convection lead to fluid flow across the den- 
dritic growth front, i. e . , normal to the heat flow direction. 
As noted, this blurs or obliterates distinction between finger 
and plume/channel patterns, leading to relatively long-range 
mixing throughout much of the solidification time. In terms 
of modeling, this more general or uniform mixing actually 
simplifies treatment of particle dispersion and survival, 
avoiding the necessity to consider local patterns and fluc- 
tuations. In fact, to complete a feasible modeling sequence, 
it may prove to be the only recourse, i.e., it may be nec- 
essary to ignore the distinction between finger and plume 
flow, in order to obtain a viable result. There are both gains 
and losses in this exercise. 

Finally, it is appropriate to remark, again, that current 
modeling procedures are inevitably substantiated by com- 
parisons with thermal analysis data by matching latent heat 
evolution to experimental cooling curves. The point to be 
reiterated is that such correlations are not a very sensitive 
function of the assumed or actual nucleation density, so that 
prediction of grain size by this method cannot be precise. 
This is the reason why some attempt should be made to 
physically describe where the nuclei originate. 
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